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The interplay between the superradiant emission of a cloud of cold two-level atoms and the radiation pressure 
force is discussed. Using a microscopic model of coupled atomic dipoles driven by an external laser, the radiation 
field and the average radiation pressure force are derived. A relation between the far-field scattered intensity 
and the force is derived, using the optical theorem. Finally, the scaling of the sample scattering cross section 
with the parameters of the system is studied. 
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1. Introduction 

When a single atom is illuminated by a laser, it scatters the light isotropically, up to polarization 
effects. The scattering process results in a force proportional to the number of photons scattered. 
Indeed, as an atom absorbs a photon from the laser of wave vector ko, it acquires a momentum 
/iko, but the average momentum change during the emission process is zero as it re-emits the 
photon in a random direction. 

As a collection of atoms, the picture changes drastically as the atoms cooperate to scatter 
the light. The phenomenon was first identified by Dicke ([7]) and it was shown to lead to a 
directional emission of the light, due to the synchronization of the atomic dipoles with the laser. 
The collective effects becomes even stronger as the atomic medium becomes optically dense and 
the radiation of the atoms starts to alter significantly the wave propagation. Among the other 
collective effects that arise, one can mention the collective Lamb shift (00, Mie resonances Q, 
as well as a reduction of the radiation pressure force (0 . 

Since the radiated light results from the interference of the waves emitted by each dipole, the 
simple relation between emitted photon and atomic recoil is lost. For example, a striking feature 
of cooperativity is the modification of the atomic recoil due to the presence of the neighboring 
atoms (l?^ i^, an effect that cannot be deduced from single-atom physics. 

We here discuss the particular relation between the directional superradiant emission, and 
the reduction of the radiation pressure force. The atomic cloud is described as a microscopic 
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ensemble of coupled atomic dipoles, and both the radiated field and the force are expressed as 
a function of these dipoles. The optical theorem is derived in this framework, and is shown to 
lead to a direct relation between intensity scattered and radiation pressure force for the cloud 
center-of-mass. 



2. Cooperative scattering model 



The atomic cloud is described as a system of two-level {g and e) atoms, with resonant frequency 
uja and position r^, that are driven by an uniform laser beam with electric field amplitude Eq, 
frequency ujq and wave vector ko = {ujQ/c)ez. The laser-atom interaction is described by the 
following Hamiltonian: 



H 



ha 



N 
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j=i k 
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where = cLEq/H is the Rabi frequency of the incident laser field and Aq = — oJa is the 
detuning between the laser and the atomic transition. In Eq. ([T]) dj = \gj){ej\ is the lowering 
operator for j— atom, ak is the photon annihilation operator and gk = {d'^oJk/2heoV)^^'^ is the 
single-photon Rabi frequency, where d is the electric-dipole transition matrix element and V 
is the photon volume. The special case where a single photon (mode k) can be assumed to be 
present in the system, was extensively investigated in Refs. 0, [77?|) . and later extended to 
include a low- intensity laser in Ref. (0, [771 HM- The system atoms+photons is then described 
by a state of the form {IS): 
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The first term in Eq. Q corresponds to the initial ground state without photons, the second 
term is the sum over the states where a single atom has been excited by the classical field. 
The third term corresponds to the atoms that returned to the ground state having emitted a 
photon in the mode k, whereas the last one corresponds to the presence of two excited atoms 
and one virtual photon with 'negative' energy. It is due to the counter-rotating terms in the 
Hamiltonian ([T]) and disappears when the rotating wave approximation is made. In the linear 
regime a ~ 1 and in the Markov approximation, valid if the decay time is larger than the 
photon time-of-flight through the atomic cloud, the scattering problem reduces to the following 
differential equation (f77| [T^ [7^ 



with initial condition /3j(0) = 0, for j = 1,...,A^. Here, F = Vg\k^l'KC = (i^k^jlT^e^fL is the 
single-atom spontaneous decay rate. The kernel in the last term of Eq. ([3| has a real component, 
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— (r'/^)'^rn=/=j[sin{xjm)/xjm] (where Xjm = ko\rj — rm\), describing the collective atomic de- 
cay, and an imaginary component, i{T/2) '^rn=/=j['^^^i^j'm) / ^jm], describing the collective Lamb 
shift ( |j-^[4-?6p . Notice that while Eq. ^ is here deduced from a quantum mechanical model, it 
can also be obtained classically, treating the two-level atoms as weakly excited classical harmonic 
oscillators ([71 [7^. 



3. Radiated field 



The radiation field operator evolves according to the following Heisenberg equation 

N 



da\f_ 1 
dt ih 



m=l 



-ikr„ 



(4) 



where the fast oscillating term proportional to exp[i(ti;fc + uja)i\ has been neglected. The scat- 
tered field is obtained by performing the sum over all the modes, considering only the positive- 
frequency part of the electric field operator 



(5) 



where £k = (/ia;fc/2eol^)^/^. Integrating Eq. (|4|) with respect to time, with ak(0) = 0, inserting 
it in Eq. ([5]), and assuming the usual Markov approximation, one obtains (j 
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When applied on the state ([2]), neglecting virtual transitions, it yields £'^1^') 
£'sexp(-iwot)|5i ...gN), with 



gifco|r-r™| 

m=l 



(7) 



Hence, the radiated field appears as a sum of spherical waves radiated by the atomic dipoles. In 
the far-field limit, one has A:o|r — r^l ~ kor — k • r^, with k = /co(r/r), so the field ^ radiated 
in a direction k reads 



fcp „ikor 
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The emitted intensity in a direction k is then derived as 
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Integrating this intensity over all directions leads to a radiated power 

- ^ 2. Wm\ + 2^ P,P^ _ ) , (11) 

\m=l rriyt^j ■' 



where we have used the equality 



kf)\d\ 



In Eq. (11), the first term corresponds to the incoherent sum of the single atom radiated power. 
The second term is an interference term; in the limit of a cloud small compared to the wavelength, 
the dipole moments have the same phase and this latter term is responsible for a superradiant 
build-up of the radiated power oc A''^ (see, e.g., Ref. ([7)). 



4. The radiation pressure force 

As for the radiation force operator acting on the jth atom, it is derived from Eq. ([T| as 

Fj = -Vr,F = Fai+Fej. (13) 

A first contribution associated to the absorption of photons of the pump appears (0 

F,, = i/iko^ {a,e*(^°*-'^-'-^) - h.c.} , (14) 
whereas the second contribution comes from the emission of the photons in any direction k: 



In Eq. (15), the counter-rotating terms proportional to exp[±.i{u)k + u)a)t\ were neglected. 

As we are interested in comparing the radiation pressure force to the single- atom case, we 
define the average radiation force F = (l/A^)X]j-^i ~ {Ftot/N)^z that measures acceleration 
of the cloud center-of-mass given by slcm = F/m, with m the single-atom mass. Note that 
this average force is N times smaller than the total force Ftot acting on the whole cloud of 
atoms. Since we consider clouds with rotational symmetry around the laser axis, this force is in 
the same direction as the incident field wave vector kg = kQe^- This average force is measured 
by time-of flight techniques in cold atomic clouds released, for instance, from magneto-optical 
traps (MOTs) and has recently revealed cooperative effects in the scattering by extended atomic 
samples (0, [73). As the scattered radiation, this force is an observable that contains signatures 
of the cooperative scattering by the atoms (0 [^. The average absorption force along the z-axis. 



resulting from the recoil received upon absorption of a photon from the incident laser, reads 

N 



(16) 



Similarly, the average emission force writes Fg = (1/A^) -^ej- Inserting the expression for at 
from Eq. Q into Eq. (15), and approximating the discrete sum over the modes k by an integral. 
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Figure 1. Scattered intensity \f{d,<f>)\'^ as given by Eq. ( |22| and Eq. | |23|l for a spherical cloud with a Gaussian density 
distribution of root mean square 10/fco- The number of scatterers is N = 1000, the detuning Aq = 0, and the incident field 
is coming from the left. The color-coded intensity is represented in log-scale. One can clearly see in red the strong forward 
emission of the sample, reminiscent of Mie scattering by clouds large compared to the wavelength. In the other directions, the 
scattered field is speckle-like due to the randomly positioned two-level scatterers, and describes the spontaneous emission of 
the cloud. Performing configuration averages would smooth out these fiuctuations, except in the backward direction where, 
in the multiple scattering regime, the well known coherent backscattering cone is recovered ^211 \22t . 



it is possible to obtain, as it was done for the radiation field operator Eg of Eq. ([7|), the following 
expression for the average emission force along the z-axis (0): 
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Neglecting virtual photon contributions, the expectation values of the absorption and emission 
forces for state ^ are 
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where we used the identity 
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ji(z) here refers the first order spherical Bessel function. Note that the decomposition into 



absorption (18) and emission (19) forces is fully compatible with classical expressions of the 
optical force (|Tgj) . where the force arises as the product between the atomic dipole and the total 
field (flO|l (i.e., including the radiation from the other atoms). 
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5. Optical Theorem 



Let us now discuss the formulation of the optical theorem in the framework of collective scat- 
tering. In the far-field limit, the field is 



gifc^^^^far)(j^^ 
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where the structure factor for the scattered field / is given by 



(21) 
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(22) 



As a consequence, the scattered intensity at a large distance r from the cloud is 



Is = lo- 



while the total scattering cross section reads 



a 



■5™ — r 2 , 



2n f'TT 
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(23) 
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Figure [T] shows the emission diagram of the scattered field and for resonant excitation and a 
spherical cloud of atoms with a Gaussian density distribution. The energy conservation imposes 
that 



(^ext — ^sca ~l~ ^abs 



(25) 



where aext and aabs are the cross sections for extinction and absorption, respectively. The ex- 
tinction cross section is then obtained from the optical theorem. In the forward direction {9 = 0) 
the total field is 



Efwd 
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(26) 



Hence, observing the field in a plane far from the atoms, the radius expands as r ~ z+^x'^+y'^) /2z 
and one obtains 
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so the intensity reads 
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(28) 



The measured intensity is the incident intensity minus the extinction intensity; in Eq. (28), the 
integration over x, y gives a factor linz/ko, and one gets 
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O-ext = -V2l™[/(0)] 



(29) 
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Hence, from Eq. ( 25 ) one obtains the relation 

f27r 
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(30) 



In our microscopic description of the light-atom interaction there is no absorption, so that 
f^abs = 0. An illustration of the validity of the optical theorem is given in Figure [2] for resonant 
light scattering by a slab containing two-level scatterers with a uniform density distribution. 



From Eqs. (22) and (29), we obtain the relation 
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Consequently, using Eqs. (18) and (19), the average force along the z-axis reads: 
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(32) 



We observe from Eq. (32) that the average radiation pressure force is not merely proportional 



to the excitation probability, i.e. l/^jPi but it is the result of an interference between the 
different atomic dipoles /3j. For this reason a measurement of the force captures the coherence 
properties of the scattering process as well as the detection of the light intensity. To make this 
point more explicit, using Eq. ([9]), it is possible to write the force as 
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(33) 



where the scattered far-field intensity is Is{6,(j)) = 2ceo|-E<j(0, (^)p. This highhghts the fact that 
the radiation pressure force, that pushes the atoms along the direction of the incident beam, is 
proportional to the net radiation flux of the scattered intensity. 

In the case of an isotropic emission (e.g., single-atom case, or cloud much smaller than the wave- 
length), the scattered intensity Ig is independent on the direction and we get Fz = {A-kt"^ /{Nc))Is'- 
the direct proportionality between scattered power and radiation pressure force is recovered. The 
cooperative effect of light scattering in such small samples is then encoded in the total scattered 
intensity Is- In the case of superradiant scattering for larger samples, a pronounced emission 
into the forward direction decreases the radiation force, as observed for example in Ref. (f^. 



6. Scaling of the scattering cross section 

In this section we are interested in understanding how the scattering cross section scales with the 
parameters of the system. We consider the case of a slab with uniform density distribution. The 
slab contains atoms and its size along the x, y, z axes is denoted by L^-, Ly, Lz respectively. 
The numerical simulations presented in figure [3] show how the scattering cross section depends 
on the optical thickness of the cloud bo = QttN / {k^LxLy). For dilute clouds of atoms we find: 



2.15 X L,j.L 



exp 



0.63 ho 
2.15 



(34) 



When the slab is optically thick, i.e. bo ^ 1^ we observe that the cross section appears to 
approach 2 x L^Ly. This factor of two corresponds to the well-known "extinction paradox" 
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Figure 2. Illustration of the optical theorem. Left: the scattered intensity integrated along (p, i.e., g{9) = J^^ dif> \ f{9, <t>)\'^, 
is shown for resonant light Aq = and a slab geometry with a uniform density distribution. The number of atomic 
scatterers is varied between 1 and 5000 (from inside to outside curves). The transverse size of the slab is Lx,y = 80/fco 
and the longitudinal size is varied such that Lz = {20 /ko)N/ 5000. This procedure allows us to vary the optical thickness 
bo = QttN / {k'^L^Ly) between 3.10~^ and 15 while maintaining the atomic density constant. The incident field is coming 
from the left and the intensity is plotted in log-scale. In addition to the forward Mie-like lobe, a lobe is also observed in the 
backward direction which we attribute to light reflection due to the sharp variation of optical index when the light hits the 
slab. Right: the blue circles represents the total scattering cross section obtained by integrating the emission diagram over 
and (j>, i.e., Osaa = I/'sq x Jq dd sin(9)g{9). In our microscopic model, there is no absorption so that (Jabs = 0, leading to 
o'ext = fsca- The optical theorem Eq. | |29| l can thus be written as o-gca = ~(47r/fcQ)Im[/{0)], which is plotted in magenta. 
The good agreement between the two curves illustrates the validity of the optical theorem. 




Optical thickness Opticai thickness 



Figure 3. Scaling of the scattering cross section. Left plot: following the same procedure as the one described in figure [2] 
we compute the scattering cross sections for different slab geometries. The results are shown in scatter plot with different 
colors. The parameters of the simulations are reported in the legend of the figure. By fitting the data, constraining the slope 
in the limit bo — > ( right plot), we obtain a scattering cross section that scales with the optical thickness bo of the slab 
according to Eq. |34| (magenta full line). 



(1^31 \24\ l for which the extinction cross section is twice as large as the one predicted by geometrical 
optics due to the diffraction contribution. The residual deviations from the factor of 2 between 
the scattering and geometrical cross sections might be associated to a still moderate size of 
our sample. For spherical dielectric spheres, aext shows an oscillatory behavior around 2ageo 
{o'geo = LxLy for our square geometry), which is damped for increasing sizes of the sphere 
l\25[ \26]i . Furthermore, in the limiting case bo <^ 1 the scattering cross section can be written 
as asca = 0.63 x {LxLy)bo = 0.63 x Nao, where ctq = 3A^/(27r) is the resonant scattering cross 
section for a single atom, whereas for off-resonant excitation (Aq = 3r), we have found that 
asca = LxLyb^Ao) (where 6(Ao) = &o/(l + ^Aq/F^)). A similar deviation in the case of a finite 
number of resonant dipoles has been noted in ( [^7| ). 
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7. Conclusion 



Wc here discussed the superradiant emission of a cloud of cold atoms, when the interference 
of the waves radiated by the atomic dipoles builds up a coherent emission. Despite the fact 
that the simple relation between absorbed photons and radiation pressure force existing in the 
single-atom case was lost, the optical theorem allowed to recover a simple relation between the 
total scattered intensity and the displacement of the cloud center-of-mass. The measure of the 
force of the center of mass of the atomic cloud contains (partial) information on the scattered 
intensity, even for large values of optical thickness of the cloud. We have computed the total 
scattering cross section which approaches a value close to twice the geometrical cross section 
of the sample, in line with the well-know extinction paradox. Further studies will include the 
precise role of the laser-atom detuning Aq, as our results indicate that 6(Ao) does not seem to 
be an universal single scaling parameter. 
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